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We calculate the photoionization with excitation- to photoionization ratios R n e and R n = Y,eR n e 
for atomic helium and positive heliumlike ions at intermediate values of the photon energies. The 
final state interactions between the electrons are included in the lowest order of their Sommerfeld 
^— i , parameter. This enables us, in contrast to purely numerical calculations, to investigate the roles 

• of various mechanisms contributing beyond the high-energy limit. The system of the two bound 

electron is described by the functions obtained by the Correlation Function Hyperspherical Harmonic 
Method. For the case of heliumlike ions we present the high energy limits as power expansion in 
inverse charge of the nucleus. We analyze the contribution of excitation of states with nonzero 
i-C ' orbital momenta to the ratios R n . In the case of helium our results for R„ are in good agreement 

with those of experiments and of previous calculations. 

^ 1 I. INTRODUCTION 

cr 

Atomic photoionization accompanied by excitation and double photoionization are much studied in experi- 
mental and theoretical works in connection with many-electron problem. In the case of two-electron systems 
■ which is considered in the present paper, the three-body Coulomb problem is the subject of investigation. 
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Most attention in experimental studies was focused on the process of double photoionization. Photoionization 
with additional excitation have not been investigated in detail. Energy dependence of the cross section ratios in 
a broad interval of the photon energies u>, dependence of these ratios on the value of nuclear charge Z, branching 
ratios for excitations of n£ subshells of a shell with the principle quantum number n are still the subjects of 
future experiments. As it stands now, there are experimental data only for atomic helium. The intensity of 
excitation of n-th shell relative to the main photoline n = 1 was measured in [1] for the photon energies up to 
several hundreds eV for the values of n < 6. A few measurements of 2 s and 2p excitations at smaller values of 
the photon energies have been carried out earlier - see [1] for references. 
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Theoretical investigation of the process requires the knowledge of the wave functions describing two electrons 
in the field of the nucleus. In initial state both electrons are bound by the nucleus. In the final state one of 
the electrons belongs to continuum, while the second one is in excited bound state. Certain approximations 
(models) for the wave functions are required. Somewhat different approximations are reasonable in different 
regions of the photon energy u>. 

We use the terminology which is similar to the one employed for much studied double photoionization [2] . It 
is known that the ratios 

cr+(w) 

of the cross sections <j£*(u>) for ionization with excitation of the second electron to the n-th level, to those 
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without excitation o\(ui) do not depend on the photon energy in high energy limit [3, 4] 

R n (oj) = const (2) 
for u — > oo. This requires anyway that uu exceeds strongly the values of single-particle ionization potentials I 

w > I . (3) 



Analysis of [3, 4] have been carried out by employing the nonrelativistic functions for description of the 
outgoing electrons. It was shown in [5] that asymptotics of the ratio R(u>) remains the same in the whole region 
(3) including the photon energies corresponding to relativistic outgoing electrons. Recall that this is not true 
for the double-to-singlc photoionization ratios [6] . 

By high energies we mean that part of the region (3) where the ratios exhibit behavior described by Eq. (2). 
At low energies the ratio I/uj cannot be treated as a small parameter. By intermediate energies we mean the 
values of the photon energies, where inequality (3) is true, while deviations of the cross section ratios from the 
high energy limit are noticeable (with the relative deviations exceeding 10 percent). For atomic helium this is 
the region from 300-400 eV till 2 keV. In the systems bound by the nucleus with the charge Z the limits of the 
interval are proportional to Z 2 . 

Since the ionization with excitation is a three-body problem, certain approximated wave functions for both 
initial and final states are required. It was shown in [3] and [4] that in the high energy limit the final state 
interactions (FSI) between the electrons can be neglected. This simplifies the problem of the description of the 
final state (under a proper choice of gauge of gauge interactions of the outgoing electron with the nucleus can 
be neglected as well). In [3] the high energy limit of the process was expressed in terms of the initial state wave 
function ^ i(r\,V2)- The high energy limits of the ratios R n for atomic helium were calculated in [7, 8] and 
dependence on the choice of the approximate function was traced. The calculations of [6] include also the Z 
dependence of the high energy limits of R n . Results of the high energy calculations for Li + are presented in [9] 
and [10]. 

At low energies there is no small parameter. All the interactions involved should be treated as accurately as 
possible. In this energy region one must make a choice of both initial and final state wave functions. The low 
energy calculations of the cross sections a^*(tu) have been carried out in [14, 15] for He and in [14, 16] for Li + . 
The paper [16] contains also results for partial cross sections <r+/(w) of ionization accompanied by excitation of 
the remaining bound electron to the subshclls with quantum numbers n and I . Low energy calculations for the 
two-electron ions with larger values of Z were carried out in [14]. 

In the papers [13, 14] the intermediate energy region was approached by extension of the low energy calcula- 
tions to this energy interval. In the present paper we move from the high energy region by including next to 
leading order of expansion in powers of u' 1 . This is achieved by inclusion of the interaction between the final 
state electrons in the lowest order of perturbation theory 

We find several attractive points in such approach. It provides the possibility to clarify the role of various 
mechanisms (in a fixed form of electron-photon interactions) representing their contributions in terms of certain 
characteristics of the initial wave function. Within the framework of the approach one can estimate the magni- 
tude of the neglected terms, thus controlling the accuracy. At the lower limit of the intermediate energy region 
numerical calculations with certain models for the final state wave functions are more precise. The discrep- 
ancy between the results obtained in numerical and perturbative approaches should diminish as u> increases. 
Hence these two approaches should supplement each other. Similar analysis of the intermediate energy double 
photoionization have been carried out earlier [17]. 

We expect the approach developed in the present paper to be useful also because of certain discrepancies 
between experimental data for helium [1] and theoretical results. The high energy limit of the ratio R n ex- 
trapolated from the data obtained in [1] is in perfect agreement with the calculated one only for n = 2. The 
disagreement increases with n reaching a factor of about 2 for n = 5. There is also visible disagreement between 
theoretical and experimental results for R n {uu) at uu <~ 200 — 400cV for n = 2,4, 5 [13, 14]. Note also that in 
the case of helium there is a discrepancy between the calculations employing various approaches (see, e.g. [13] 
and [14]). Moving from the high energy region can be instructive also since (as noted in [14]) the R matrix 
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approach, widely used in the low energy calculations becomes unstable at the high energies. Finally studies of 
Z dependence of the ratios R n may be of interest in connection with increasing attention devoted to physics of 
the multicharged ions. 

We calculate ratios (1) of photoionization accompanied by excitation of the residual ion for helium atom and 
light hcliumlike positive ions. We obtain also more detailed characteristics 

Such ratios are also detected in the low energy experiments [18]. The ratios defined by Eq. (1) can be represented 
as 



In this paper the calculations are carried out with inclusion of next-to-leading order terms of expansion of 
the ratios (1) in powers of iv^ 1 . This means that for ns states we calculate the high energy limits of the ratios 
(1) and the correction of the order For nt states with / >lwc obtain the leading order of expansion in 
l/ui. 

In the limit (3) all the interactions of the outgoing electron can be treated perturbatively [3]. In the high 
energy limit of final state interactions (FSI) of the outgoing electron with the electron bound in the residual 
ion can be neglected. The excitation following photoionization is due to the specific correlation in initial state 
known as shake-up (SU). Only s states can be excited by this mechanism. Excitations of the states with nonzero 
values of angular momentum i are quenched by a small factor of the order I/oj. 

We can present the ratios (4) as 

R ns (u) = A n + ^B n0 (5) 
(with I the electron binding energy in hydrogen) for I = 0, while for I > 1 

Rne(ui) = —B ne , (6) 

In the atomic system of units used through the paper (e = m = h = 1, c = 137) Io = 1/2. 

The coefficients A n and Bg n with I > 1 do not depend on the photon energy while we show that B n0 contains 
a smooth dependence on lo. Now we can present the ratios R n defined by Eq. (1) as 

R n (uj) = A n + —B n ; B n = B n j . (7) 



In the SU mechanism the interactions of the outgoing electron with nucleus can be treated perturbatively [3] . 
Excited electrons can be described by the Coulomb field wave functions. Thus, all the specifics of this three- 
body problem is contained in the wave function of the initial state. The ionized electron approaches the nucleus 
at the distances which are much smaller than the size of the atom. The SU cross section is thus determined by 
initial state wave function ^E , i(r 1 ,r 2 ) at electron-nucleus coalescence point, i.e. by ^(ri = 0,r 2 ). 

The SU probabilities depend on n in terms of the wave function and of momentum p n of the outgoing electron 

pl = 2e„ 

with e n being the energy of the outgoing electron. In the lowest order of expansion in powers of I/u) we can put 



2 2 2 o 

Pn = Pl = V =2(J. 



(8) 
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The SU terms with this value of p determine the high energy limit of the ratios 

lim R nl {uj) = lim R ns (uj)S m = A n . (9) 

Now we consider three types of contribution beyond the high energy limit, in the same way as it was done in 
[17] for the double photoionization. The kinematical corrections to SU ratios are caused by taking into account 
n dependence of momentum p n in the SU amplitudes. This provides the contributions to the terms B n on 
the RHS of Eq. (5). Note that these corrections are proportional to the small parameter I/oj, containing also 
dependence on specific parameter 7r£ z with 

iz = - ■ (10) 
P 

One has to have in mind that corrections of the order 7r£z/w drop as uj~ 3 / 2 but contain a numerically large 
coefficient. We shall not treat 7r£z as a small parameter, but include it exactly. The dependence of the cross 
sections on n^z known to be presented by the Stobbe factor S(n^) = exp (— 7r£) which is common for the 
photoionization processes [19, 20]. These corrections are expressed in terms of SU contributions A n to the 
ratios (5), which appear only in the ratios R ns . 

In the next to leading order the excitation energy can be transferred to the second electron also by the initial 
state interactions (ISI) beyond the SU. In this case the terms of the order l/u> and Z 2 /to come from the higher 
terms rf/r 2 and (rir 2 )/r 2 of the expansion of initial state function 5 , i(r 1 ,r 2 ) at n — > 0. Thus the contribution 
will be presented in terms of the derivatives of the initial state wave function, integrated with the Coulomb field 
function of the bound state. 

The excitation energy can be transferred also by the final state interactions (FSI) between the final state 
electrons. We include the FSI by perturbative method developed in [21]. The FSI amplitude is presented as 
power series of the Sommerfeld parameter of the interaction between the fast outgoing electron and that of the 
residual ion 

z = l (") 

while v is their relative velocity. Thus, the square of the amplitude is presented as power series in £ 2 = l/2e 
with e being the energy of the outgoing electron. Looking for the terms of the relative order lo^ 1 in the cross 
sections, we must include the lowest correction of the order £ 2 , putting 

« 2 = h- (12) 

The FSI contributions are presented in terms of matrix elements of relatively simple operators sandwiched by 
the function ^ / i(ri,r 2 ) and the Coulomb function of the electron in residual ion. The states with any angular 
momenta I can be excited by the FSI in the next to leading order of uj^ 1 expansion. 

Thus, all the contributions up to the order w _1 will be presented in terms of certain characteristics of 
the initial state wave functions. We employ the functions obtained by Correlation Function Hyperspherical 
Harmonic Method (CFHHM), obtained in [22]. The CFHHM functions have been employed successfully for 
investigation of the parameters of the bound two-electron systems [23] and of some characteristics of the double 
photoionization. Also, the method of inclusion the FSI [21] enabled earlier to remove the discrepancy between 
experimental and theoretical results in creation of vacancies in electronic shells during nuclear transitions and 
in single photoionization of the p states. In [17] it was used for investigation of the double photoionization at 
intermediate energies. In the present paper we use the CFHHM functions and the perturbative treatment of 
FSI for investigation of photoionization with excitation. 

Note that for the system containing larger number of electrons the picture is more complicated. Considering 
ionization with excitation of the subshell with £ = 1 we find for ionization without excitations a + <~ oj~ 9 / 2 , while 
for ionization accompanied by excitation to an s state er+* ~ uu~ 7 / 2 (the ISI provides admixture of two s state 
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electrons to the system containing two p state ones) . Hence the corresponding ratio increases proportionally to 

OJ. 

Our analysis is completely nonrelativistic. We neglect the terms of the order u)/m in the wave function of 
the final state, and in the operator of the photon-electron interaction. The latter means that we are using the 
dipole approximation. We assume also (Z/137) 2 <C 1, to neglect relativistic effects in the initial bound system. 



II. GENERAL EQUATIONS 



The cross section of photoionization accompanied by excitation of the residual ion into a state with the 
quantum numbers n, £, m can be written as 

d<; = j| E \^\ 2 ^ - *» - h) ^ • (13) 

Here Ii denotes the ionization potential of K electron in initial state atom. The factor 2 is due to two electrons 
in the K shell. The ovcrline shows that the averaging over the photon polarizations have been carried out. The 
angular dependence of the amplitudes can be written explicitly due to the dipole approximation employed The 
amplitude 

F nlm = ($„ £m | 7 |*) , (14) 

with 7 being the operator of interaction between the photon and an electron, while W and <& n im describe the 
initial and final two-electron states, can be represented as 

F nlm = {^) 1 ' 2 ^^T nlm . (15) 
c 

After averaging over the photon polarizations one obtains 

<t = |- JcEiW- (16) 

m 

If the FSI are neglected, the final state function is 

®ni m {ri,r 2 ) = i)f{p n ;r 1 )ip nlm (r 2 ) (17) 

with the functions ipf and ip n ,e,m being just the continuum and bound state single-particle wave functions in 
the Coulomb field. If condition (3) is valid, the interactions of the outgoing electron with the nucleus can be 
included perturbatively. Using the velocity gauge for the operator 7, i.e. 

7(7*) = —i(e ■ V) 

with e standing for the photon polarization, we obtain the leading contribution of expansion in powers of p~ 2 
as coming from the plane waves. Following [17] we can separate two scales in the interactions between the fast 
outgoing electron and the nucleus. Those taking place at small distances of the order «r c with r c = l/Z 
being the characteristic size of the atom are expressed in terms of the parameter tt£,^ = Z/p( n \ Such 

contributions can be calculated explicitly, producing the factor 

N^) = N r (^)e-^ (18) 

with N r {s}^) = {2ttS}^ /(l — cxp (— 27r£^))^ being the normalization factor of the nonrelativistic Coulomb 

continuum wave function. The interactions which take place at the distances of the order r ~ r c can be presented 
as p~ 2 series thus cancelling in the ratios (1) and (2). Thus we can put 

*SL(n,r a ) = M^V^Wra). (19) 
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Following [22] we present the factors N 2 {S}^ ) in the cross sections as 

N 2 (^ l) )= h(*zP)e-«& (20) 

with the function h{(}^) — 27r£^/(exp (7r£^) + ex P ( — ^z^)) containing only weak dependence on parameter 
7r£^\ Thus we can put h(n^^) — h(ir£z), with £z defined by Eq. (10). Hence, 

N 2 (4 n) ) = h{nZ z )e-*& . (21) 

The second factor on the RHS of Eq. (21) is the Stobbe factor mentioned above. 

We shall present the perturbative FSI contributions also in terms of the function (17). Thus the ratios (1) will 
be presented in terms of the matrix elements of initial state two-electron function and the Coulomb function of 
the excited electron. 



III. AMPLITUDES BEYOND THE SHAKE-UP 

Following the analysis given above, we present the amplitudes for ionization with excitation beyond the SU 

as 

-Tnim — r nl m T r n g m ~r r n lm X^^J 

(s) 

with F n g m standing for SU amplitude, which includes kincmatical corrections to the high energy limit, while 
^nlm an d Fnlm are ^ ne contributions caused by correlations in ISI and FSI correspondingly. 



A. Amplitudes without inclusion of final state interactions 

Since in this subsection we neglect the interactions between the outgoing electrons, we can use Eq. (14) for 
the amplitude with the final state wave function presented by Eq. (17). Recall that we use the operator 7 in 
the velocity form. This provides 

Frum = (4 7 r) 1 /2( e .pJ7V(^) J d 3 r 2 r nem (r2)*i( Pn ;r 2 ). (23) 

Here 

*i(p n ;r 2 ) = J d 3 r 1 ^ i (r 1 ,r 2 )e- i( -Pn r ^ (24) 

is the partial Fourier transform of the initial state wave function in variable r\. 

Since the integral over r 2 on the RHS of Eq. (23) is saturated at r 2 <~ r c , while p n » r^ 1 , we need expansion 
of the function ^ (p n ; r 2 ) in inverse powers of p 2 n . It is convenient to employ the Lippman-Schwinger equation 

%(p n ,r 2 ) = — J(p n ,r 2 ); J( Pn ,r 2 )= / — e^^r, r 2 ). (25) 

Pn J r 

The integral on the RHS is dominated by r ~ p~ x <C r 2 . Thus the expansion in p~ 2 can be carried out by 
expanding the function ^(r, r 2 ) in powers of r in the integrand on the RHS of Eq. (25). 
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1. Shake-up with kinematical corrections 

Presenting 



J( Pn ,r 2 ) = lim| ^ e -<(Pnn-*r*. (PjP2)} (26) 
we obtain for the leading order contribution 

J(p n ,r 2 ) = J*(0,r 2 ), (27) 
which enable to write for the SU amplitudes 

F nL = a{p n )S n 8 m 5 m0 ; S n = (An) 1 / 2 J dr 2 r 2 2 ^(r 2 )H>(0, r 2 ) (28) 

with 

a(p n ) = ( e p n )N(n^)^§, (29) 

P Pn 

while the upper index (r) in Eq. (28) labels the radial part of the Coulomb function ip n o- In the leading order 
we should neglect the dependence of p n on n, putting p n = p, just as in Eq. (8). The high energy limit of the 
amplitude (28) is thus 

Fill = <P)S n - (30) 

In the next to leading term we must include the n dependence of p n . Since the residual ion contains only one 
electron, the latter is described by the Coulomb wave function, and thus 

pi = p 2 - 25 n , (31) 

where 

is the excitation energy of the electron in the final state ion. 



2. Contributions of correlations in the initial state 

Now we return to Eqs.(24), (25), looking for higher order terms of expansion of the function \ti(r,r 2 ) at 
r — > 0. Since the CFHHN functions are expressed in terms of variables r — \r\,r 2 — \r 2 \,p = \r — r 2 \, we 
present the expansion in terms of the function ^(r, r 2 ,p) = v Pi(r,r 2 ) 

*(r,r 2 ,p) = (l + r i Vi + ^rir j ViV J -)*(r,r2,|r-r 2 |). (33) 

Here we put *(r, r 2 , \r — r 2 \) = ^(O, r 2 , r 2 ) after the derivatives are calculated. Using Eq. (33) we find that at 
small r 

(r, r 2 ,p) = *(0, r 2 , r 2 ) + r<(r, r 2 , r 2 ) - rr^(0, r 2 ,p) + 

r 2 n r 2 (l-T 2 ) , 
+ — * (r,r 2 ,r 2 ) + — itf (0,r 2 ,p) + 

2 T 2 

r 2 T 2 „ 

+ -^-* P (0, r 2 , p) - r 2 r* rp (r, r 2 ,p). (34) 
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Here r = (r ■ r 2 )/r ■ r 2 . The derivatives ^l' r and ty' (and those of the second order) are taken at the points 
r = and p = r 2 . The higher terms of expansion in r contribute to the higher order corrections in l/p to the 
amplitudes. Thus, they are neglected. While evaluating the next to leading order terms we must put p n = p. 

Using Eqs. (26) and (34) we find nonzero contributions to the amplitudes with the angular momenta 1 = 
and 1=1. For excitation to s states we obtain 



F (2) 
r n00 



a(p)Q„^; 



Q n = -(4^) 1 / 2 | dr 2 rl^{r 2 ) x 



K(r,r 2 ,r 2 ) + ^(0,r 2 ,p) + ^-& p (0,r 2 ,p) 



(35) 



with the upper index (r) labelling the radial part of the single-particle Coulomb field function as in Eq. (28). 
The function a{p) is determined by Eq. (29). Other notations are explained in the text below Eq. (26). For 
excitation into p states, choosing the direction of the outgoing electron momentum as the axis of quantization 
of the angular momentum, we obtain 



F$ m = ia(p)P n S -f; P n = (4*)V» ^ J dr 2 r 2 2 ^(r 2 )^(0,r 2 , p)r a . 



(36) 



Thus interactions in the initial state provide corrections of the order p~ 2 to the cross sections of excitations 
into s states. The dependence of the wave function on the interelectron distances p — \n — r 2 \, which describes 
the electron correlations enables also excitations into p states. Excitation to the states with higher orbital 
momenta due to the ground state correlations only are still impossible. 



B. Contribution of the final state interactions 



Now we include the final state interactions. Following [21, 24] we present the final state wave function as 

$W = (l + GV ee + GV ee GV ee )$W (37) 

with <&(°) being the wave function (17), where the FSI have been neglected (here we omit lower indices), V ee 
is the electron-electron interactions, G is the propagator of the system of two non-interacting electrons in the 
Coulomb field of the nucleus. The second and third terms on the RHS of Eq. (37) correspond to one and two 
interactions between the final state electrons, thus being proportional to the powers of the parameter £ defined 
byEq. (11). 

The two last terms on the RHS contain infrared divergent contributions caused by the Coulomb interactions 
V ee . It was shown in [21] that the infrared divergent terms cancel in each order of the expansion of the square 
of the amplitude |-F| 2 . The situation is similar to that with the infrared singularities in the e — N scattering 
analyzed in [25]. The cancellation can be illustrated by assuming the electron interactions to be defined as 
V ee (r) = lim^o e~ vr jr. The contributions lni^ emerge in intermediate steps but vanish in the final expression 
for |F| 2 . 

Explicit expressions which include the FSI in process with the fast outgoing electron up to the terms of the 
order £ 2 have been obtained in [21, 24]. The first order amplitude F^ x \ corresponding to the second term on 
the RHS of Eq. (37) is mostly imaginary. The real part of F^ 1 ) is suppressed by additional power of p~ x and 
thus can be written as being proportional to The second order amplitude is mostly real. Thus Im F^ 1 " 1 <~ £, 
ReF^' 1 ) - f, RcF^ 2 ) - £ 2 , ImF^ 2 ' - <f (we do not trace the dependence on Z here). The FSI amplitudes 
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can be presented as [21] 

F%2 = a(p) (i^ nlm \ ln(r 2 - r 2z )v\*i) + y <^ m |r ^H*i>) ; (38) 

= <VW| ln 2 (r 2 - r 2 ,)H^) . (39) 

Here ^ = ^(ri = 0, r 2 ), is a function of r 2 = |r 2 |, z is the direction of the momentum of the outgoing electron. 
Recall that z axis is chosen for the quantization of angular momentum, r$ is the Bohr radius. Thus all the 
contributions to F^) have nonzero values only for m = 0. 

For s states both amplitudes F* 1 and F-^ 2 are important, since the terms containing the factor £ 2 interfere 
with the SU amplitude. For s states 

F^ = a (pMU n + ev n ); F^ = -a(p)ew n , (40) 

with 

U n = (4^) 1 / 2 y"dr 2 r 2 V'i r ) (r 2 )lnr 2 z/* J (0,r 2 ); 
( 47r ) 1/2 [j 2,M, xd*i(0,r 2 ) 



Kn = 



2 

TU„ = J dr^l^ir^ir^^). (41) 

For the states with I ^ there is no interference with the SU amplitude. Thus only the first term of the 
amplitude F^ 1 ^ is important. We can present 

= iaiS nt c t + Q{e) (42) 

with 

S nt = (Air) 1 / 2 J dr 2 r 2 2 ^ e \r 2 )*(0,r 2 ), (43) 

while 



/ dr 2 r 2 ^(r 2 ) ^ r ; 



Q = ; ' ■ ' 



. I dtln(l-t)^(t) = -^ I y ) (44) 
-1 

and Vt is the Legendre polynomial. 



IV. THE RATIOS 

Now we can calculate the cross sections and the ratios (1) and (4). The cross sections are related to the 
squares of the amplitudes \F\ 2 by Eqs. (13) and (16). We start with calculation of the values of \F\ 2 . 



A. Excitation of s states 

Expressions for excitation of s states have the most complicated structure 

\F ns \ 2 = a 2 ( Pn )Sl + [2S n Q n + 2S n (V n + W n ) + U 2 ] . (45) 
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Here the first term on the RHS stands for SU contribution with account of kinematical corrections - Eqs. (28) 
and (29). The first term in the square 

brackets comes from interference between SU and ISI amplitudes - Eqs. (30) and (35). The second term in 
brackets is caused by interference of SU amplitude presented by Eq. (30) with the first and second order FSI 
amplitudes presented by Eqs. (40) and (41). The last term in brackets is a purely FSI contribution. 

In order to obtain contribution of the first term on the RHS of Eq. (45) to the ratio (1) we must include the 
n dependence of the phase volume in Eq. (13) for the cross section. As a result, for purely SU ratio we find 



R su . 



S 2 lPri 



(n) 
Z 







(46) 



with p n defined by Eq. (31). In the lowest order of expansion in powers of I/u> the dependence of R„ on 
parameter tt^z is just the same as it would result in expansion of the RHS of Eq. (46) in powers of tt£z- 
However this is not true for the higher order terms of I/w expansion. 

The contribution of the other terms to the ratios R n can be found as their ratios to squared amplitude of 
photoionization without excitation |-Fi s | 2 , where the corrections of the order I/w also should be included. This 
gives 



FL 



o2 

»M = | + - 

1S n 



i 

2lo 



^ — (1 - — )(1 

St 2 1 ni >{ 



S n 

sl 



+ u 2 

1 ^ 11. 



o2 

°n jj2 



(47) 



The first term on the RHS is the high energy limit of the ratio. Note that Eq. (47) provides exact dependence on 
parameter n^z in next to leading order of I/co expansion. Expression in square brackets on the RHS of Eq. (47) 
should be identified with parameter B n Q introduced by Eq. (5). Separating energy independent contributions 
and the terms, which depend on the photon energy through parameter tt^z we write 



B n Q = d n + TT^zfn 



(48) 



with 



fn 



Sl& 
2 



2S n 



S? 2 1 n 2j+ SI 
2 ' 1 



Qn + V n + W n 



V 1 + W 1 ))+U:-^; 



o2 

Sf 



(49) 



B. Excitation of p states 



Excitations of the states with I ^ can take place only beyond the SU approximation. The contributions of 
ISI and FSI are expressed by Eq. (36) and by the first term in brackets on the RHS of Eq. (38) correspondingly. 

F nl = i-(P n + S nl ), (50) 
P 

leading to 



(51) 
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c. 



Excitation of the states with I > 2 



In this case only the FSI contribute. The amplitude is expressed by Eq. (43) giving 



1 r 2 <? 2 



(52) 



2lo Si 



with ci defined by Eq. (44). 

Now we turn to analysis of particular cases. 



V. 



THE CASE OF HELIUM 



A. High energy limit 



In Table I we compare our results for the high energy limits with those extrapolated from the experimental 
data [1] and with the results of previous calculations [7, 9, 13, 14]. There were also several publications of the 
high energy limit for R 2s only. The pioneering calculation [26] provided i?2s = 4.61 • 1CP 2 , while the latest 
available result is R 2s = 4.79 • 1(T 2 [27]. 

One can see that various theoretical approaches provide very close values for high energy limits at n < 6. 
However the values of these limits extrapolated in [1] from the experimental data are in perfect agreement with 
the theoretical results only for n = 2. Discrepancy between theoretical and experimental results increases with 
n rapidly. 



Now we consider the contributions beyond the limit (9). Let us start with excitation of ns states. As we 
showed above, the contributions beyond the high energy limits come from kincmatical corrections to the SU 
term and from the initial and final states electron-electron interactions. One can see from Table 2 that the FSI 
and ISI contributions are positive, with the FSI term being about 3 times larger than the ISI one for all values 
of n. The kinematical corrections are negative at /i > 1, corresponding to u> < 540 cV. At larger lu values they 
become positive. Their contribution to the parameter B n o defined by Eq. (5) does not exceed 10 percent for 
oj < 1 keV. They increase with u>, becoming as large as 25-30 % in the limit \i <C 1 (w » 500 eV). In this limit 
the FSI contribution determines about one half of the parameter B n Q. Note that the relative role of the three 
contributions does not vary much with n. 

Contributions to excitation of p states come from FSI and ISI. Actually the former one dominates, providing 
more than 4/5 of the total contribution to B n \. The result of calculations are presented in Table 3. One can 
see that B n0 and B nl provide contributions of the same order of magnitude to the energy dependent part of 
the ratio i?„(w) defined by Eqs. (1) and (7). 

Thus the energy dependent parts of the ratios R n {ui) determined by parameters B n , are dominated by 
contributions of s and p states. The d states provide corrections of about 10%, while contributions of states 
with larger values of orbital momenta are negligibly small. The coefficients B n are dominated by FSI which 
provide more than 70% of the values. 

Excitation of the states with I > 2 are due to FSI only. The cross sections of d states excitations are several 
times smaller than those of p states still providing noticeable contributions to parameter B n (Eq. (71)) for 
n > 3. The relative role of excitation of d states slowly increases with increasing n - see Table 3. Excitation of 
states with higher values of t drops rapidly with increasing of I. For example, the cross section for excitation 



B. Beyond the high energy limit 
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of 4/ state is about twenty times smaller than that of Ad state. The values of B n i and B n are presented in 
Table 3. 

Energy dependence of the relative role of excitations to the states with t > 0, expressed by the ratio cr^* /er+* = 
Rnt/Rn is shown in Fig. 1. 

The ratios B n /A n converge to certain limiting values while n increases. This is due to similar n~ 3 behavior 
of both characteristics at large n. 



C. Comparison with earlier results 

Now we compare our result with experimental and theoretical data obtained by the others. In Fig. 2 we 
present the ratios R n for 2 < n < 6, calculated in the present work and measured in Ref. [1]. Wc show also the 
results of calculations carried out in [14], where the intermediate energy region was reached by moving from the 
low energies. 

We see that for n = 2 our results are in good agreement with those of [1, 14] for uj > 400 eV. As expected, 
there are noticeable deviations from results of [1, 14] at smaller values of u>. We do not show the results of 
calculations for lo ~ 900 cV obtained in [11, 27], which are also in good agreement with those of the present 
paper. For n = 3 we find a good agreement with experimental and theoretical results at all values of lo. For the 
cases n — 4 and n = 5 our results are close to those of [14], with both sets of the calculated values exceeding the 
experimental data at w ~ 300 — 400 eV. For n = 6 the experimental results are available only for u> < 160 eV, 
where the accuracy of our approach is poor since £| > 0.34. However the deviations between our results and 
experimental data is are not large for n = 6, as well as for the other values of n in this energy region. 



VI. Z DEPENDENCE 
A. High energy limit 

It was shown in [6] that the SU ratios A n drop as Z~ 2 at Z ^> 1. The tendency is illustrated by the results 
of calculations for Z < 10 presented in Tabic 4. One can see that the convergence to Z~ 2 law becomes better 
for larger values of n. In [28] Z dependence of high energy limits for double-to-single photoionization ratio 
was traced and presented as a Z^ 1 series. We can write similar presentation for the ionization followed by 
excitation. Assuming that A n can be approximated by two terms of the series one has 

^ = ^ + ^- (53) 

The values of a n and b n are given in Table 5. The convergence of Z^ 1 series is faster than in the case of double 
photoionization [28]. Also, in contrast to the double photoionization case, the leading Z~ 2 terms underestimate 
the values of the ratios. Note that in approach of [6], where all the interactions between the electrons where 
treated perturbatively, 

A n = % (54) 

with c„ also presented in Tabic 5 (a numerical error was corrected in [10]). As expected, our values of a n are 
close to c„. 

The results presented in Table 5 illustrate also the tendency to n~ 3 behavior. The values of the product 
n 3 ■ i?l s } for n = 5 and n = 6 differ by 5% for Z = 2 and by 4% for Z = 10. 
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B. Beyond the high energy limits 

In this case Eq. (3) can be written as 

Z 2 

u » — . 

In order to trace Z dependence of the characteristics, we consider the limit Z> 1. Let us start with excitations 
of s-statcs. One can see that in Eqs. (48) and (49) the ratio S\/S n <~ 1/Z, while Q n /S\ <~ Z, V n /Z <~ Z, 
W n /Si <~ 1, U n /S\ <~ 1. Thus in Eq. (5) A„ <~ 1/Z 2 , _B„ S <~ 1. Hence, we can present 

■RnO = ^1 + ^~ r ™<^ , (55) 

with r„ = £> n o/Z 2 ^4„. Using Eqs. (48) and (49) we can write 

r n o = r% + ir&rlo (56) 

with r^ = d n /Z 2 A n and 

/ fn n 2 -1 
r - 

In similar way we can present for I > 1 



r "°"^X" ~2n^- (57) 



Z 2 

Rnt = -^r nt A n , (58) 

with r„£ = B n i/Z 2 A ni while parameters B n e are introduced by Eq. (6). 
For the cross section ratios R n defined by Eq. (1) we present 

Rn = (l + |- r„) A n ; r„=r^+7r^; r^=7^ + ^rw; r£ = r£ . (59) 

For illustration we present characteristics of the process for Z = 10. Interplay of the three types of contribu- 
tions to the parameter B n0 describing excitation of s states is shown in Table 6. As well as in the case Z = 2 
the FSI provides the main contribution. However the domination is less pronounced than in the case of helium. 

As one can see from Eq. (57), the ratios r^ do not depend on Z. Dependence of parameters r n e and r n on 
Z is illustrated by the results presented in Table 7. The values of for Z = 10 are somewhat larger than for 
Z = 2. This is mainly due to larger contribution for excitation of p states in the case Z = 10. On the other hand 
the role of excitation of d states becomes smaller - see Fig. 3. The ratio of for Z = 2 and Z = 10 exhibits 
very weak dependence on n. The values of r n £ and r n converge to certain limiting values while n increases. 
This is due to similar n~ 3 behavior of the parameters B n £ and A n at large n. 

To estimate the limiting behavior of the ratios R n for Z 3> 1, note that the first term on the RHS of Eq. (7) 
for R n depends on the nuclear charge as Z~ 2 . Since the values of are several times smaller than they can 
be neglected for n^z < 1. At these energies the second term contains only weak dependence on Z. 



VII. SUMMARY 



We have considered photoionization accompanied by excitation for helium atoms and positive two-electron 
ions. We focused on the case of intermediate photon energies, for which expansion of the amplitudes in powers 
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of lo^ 1 is possible, while account of the lowest term only is not sufficient. We included the final state interactions 
between the electrons in the lowest order of their Sommcrfcld parameter. This enabled us to analyze the role 
of various mechanisms of transferring the excitation energy. 

We calculated the ratios R n e of the cross sections <r^* for ionization, accompanied by transition of the second 
electron to the bound state with quantum numbers n, I to the cross section for ionization without excitations 
a J,, and also found the sums R n = ^2 e R n e - Eq. (41). 

Following [17], we separated three types of contributions beyond the high-energy limit. These are the kine- 
matical correction to the shake-up (SU) terms, and the contributions describing the transfer of excitation energy 
by initial state and final state interactions (ISI and FSI correspondingly). The FSI were included by the pertur- 
bative approach developed in [21] and employed in [24]. This enabled us to extract the energy dependent factors, 
presenting the amplitudes in terms of the matrix elements containing the initial state wave functions. The latter 
have been obtained in [22] by Correlation Function Hypcrspherical Harmonic Method. These functions were 
employed in atomic physics calculations earlier [23]. 

We carried out the calculations, taking into account the next-to- leading terms of expansion in powers of 
The kinematical corrections to the SU terms depend also on the parameter ir£z- Dependence on this parameter 
was included exactly. 

The cross sections for excitation of ns states have the most complicated structure. In this case we had to 
include kinematical corrections to SU terms. The ISI amplitudes are proportional to and we included their 
interference with the SU amplitudes. The first and second order FSI amplitudes contain the factors i£ and 
£ 2 . Thus we had to include the interference between SU and FSI amplitudes and a purely FSI term. All the 
corrections should be included in the expressions for ionization cross sections without excitation a^ as well. 
The cross section for excitation of p states was determined by ISI and FSI mechanism, with both amplitudes 
being proportional to i£. Ionization accompanied by excitations to the state with £ > 2 took place only due to 
the FSI. 

For the case of helium we found the values of the high energy limits for n < 6 to be in agreement with 
those, calculated by the others - Table 1. For excitations of ns states we found the FSI to provide the largest 
contributions. Excitations of np and ns states, provide the contributions of the same order of magnitude to the 
energy dependent parts of the ratios R n . Excitations of nd states determine a corrections of about 10% to R n . 
Excitation of states with I > 3 are negligibly small - Tables II and III. For the atomic helium we carried out 
detailed comparison with earlier experimental and theoretical results. We found good agreement at ui > 400 cV 
for n — 2 and even at smaller values of ui for the larger values of n - Fig. 2. 

For larger Z we found approximate formula (53), which presents the high energy limits in the form of Z^ 1 
series, with the leading terms of expansion being consistent with the earlier results [6, 10]. Excitation of ns 
states beyond the high energy limit is still dominated by the FSI. The role of excitation of the states with I = 1 
increases, e. g. for Z = 10 transitions to np states provide the largest contribution to the energy dependent 
part of R n . The role of excitation of nd states drops with Z. These results are illustrated by Table VII. In the 
limit Z» 1 the ratios R n can be presented as the sums of two terms. The high energy limit term does not 
depend on ui, dropping with Z as Z~ 2 . The second term drops as varying with Z slowly. 

For the case of helium, as well as for the ions with larger values of Z, the contribution of ISI to the ratios R n 
is about 10%. Hence, the ratios R n are determined by the kinematical corrections to SU and by the FSI. 
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Figure captions 



Fig.l. Energy dependence of the relative role of excitations to the states with £ > 0, expressed by the ratio 
X n (. = cr^e / a n* = Rni/Rn for the case of helium [Z = 2). 

Fig. 2. Energy dependence of the ratios R n in 1CP 2 units for the case of helium. The dots stand for the 
experimental data of [1] . The solid lines show the results of the present work. The dashed lines show the results 
of the calculations carried out in [14]. 

Fig. 3. Energy dependence of the relative role of excitations to the states with I > 0, expressed by the ratio 
X nt = cr+,7<7+* = R n i>/R n for the case Z = 10. 
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TABLE I: The high energy limits of the ratios R n3 in percent for helium atom. The column "Theory-I" [7] stands for 
early high energy calculation with a variational initial state function, "Theory-II" [9] results obtained using many-body 
perturbation theory, Theory-Ill [13] and Theory-IV [14] are the extensions of the low energy results with multiconfig- 
uration Hartree-Fock and variational ground state wave functions correspondingly. The last column shows the result 
obtained in [1] by extrapolation of their experimental data of [1] with statistical errors given in parentheses 



n 


Theory-I 


Theory-II 


Theory-Ill 


Theory-IV 


This work 


Experiment 


2 


5.34 


4.79 


4.78 


4.79 


4.80 


4.80 (13) 


3 


0.66 


0.592 


0.605 


0.596 


0.590 


0.543 (33) 


4 


0.21 


0.19 


0.200 


0.197 


0.195 


0.118 (37) 


5 


0.100 


0.09 


0.092 


0.091 


0.0900 


0.048 (30) 


6 


0.055 


0.05 


0.050 


0.0515 


0.0493 





TABLE II: Contributions of various mechanisms to the value B n o defined by Eq. (5) for the case of helium. Here /i = ir£z; 
(j, = 1.04 or u) = 500 eV, (j, = 0.73 for u) = 1 keV. 



n 


Kinematical corrections 


ISI 


FSI 


B„0 


2 


0.072(1 - n) 


0.027 


0.094 


0.193 - 0.072/i 


3 


0.105(1 - /i)(-l) 


0.50(-2) 


1.46(-2) 


(3.01 - 1.05/t)(-2) 


4 


0.036(1 -m)(-1) 


1.76(-3) 


5.00(-3) 


(1.04-0.36/t)(-2) 


5 


0.174(1 - /i)(-2) 


0.85(-3) 


2.40(-3) 


(4.99- 1.74/t)(-3) 


6 


0.095(1 -m)(-2) 


0.47(-3) 


1.35(-3) 


(2.77-0.95/t)(-3) 



TABLE III: Parameters B n e and B n of the energy dependent contributions to the ratios R„t and R n defined by Eqs. 
(5)-(7). 



States 


B n e 


B„ 


2s 


0.193 - 0.072/x 








2p 


0.130 




0.323 


- 0.072/i 


3s 


(3.01 - 1.05/tt)(- 


-2) 








3p 


1.86(-2) 




(5.13- 


1.05/i)(- 


-2) 


3d 


3.07(-3) 










4s 


(1.04-0.36/tt)(- 


-2) 








4p 


0.62(-2) 




(1.80- 


0.36/ti)(- 


-2) 


4d 


1.35(-3) 










4/ 


3.9(-5) 










5s 


(4.99- 1.74/i)(- 


•3) 








5p 


2.93(-3) 




(0.86 - 


0.17/i)(- 


-2) 


5d 


0.70(-3) 










6s 


(2.77 -0.95/i) (- 


•3) 








6p 


1.63(-3) 




(4.80 - 


0.95/i) (- 


-3) 


6d 


0.40(-4) 
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TABLE IV: The values A n Z 2 ■ 10 2 , with A n defined by Eq. (9). 



n 


Z = 2 


Z = 3 


Z = 4 


Z = 6 


Z = 10 


2 


19.1 


14.9 


12.8 


11.4 


10.4 


3 


2.36 


2.18 


2.06 


1.93 


1.84 


4 


0.781 


0.749 


0.722 


0.692 


0.660 


5 


0.360 


0.351 


0.340 


0.327 


0.316 


6 


0.197 


0.193 


0.188 


0.182 


0.176 



TABLE V: The values of coefficients a n and b n in Eq. (53) and of coefficients c„ in Eq. (54). 



n 


a„(-2) 


M-2) 


c n (-2) 


2 


8.9 


15.0 


9.2 


3 


1.7 


1.4 


1.7 


4 


0.61 


0.48 


0.64 


5 


0.30 


0.16 


0.30 


6 


0.17 


0.08 


0.17 



TABLE VI: Contributions of various mechanisms to the value B n defined by Eq. (5) for the case Z = 10, /i = n£z- 



n 


Kinematical corrections 


ISI 


FSI 


B n Q 


2 


0.390(1 




-1) 


1.91(- 


-2) 


4.96(- 


-2) 


0.104 


- 0.039^ 


3 


0.820(1 


-/*)(- 


-2) 


0.41(- 


-3) 


1.00(- 


-2) 


(2.23- 


0.82 / u)(-2) 


4 


0.326(1 


-m)(- 


-2) 


1.56(- 


-3) 


3.85(- 


•3) 


(8.67- 


3.26/i)(-3) 


5 


0.152(1 


-m)(- 


-2) 


0.76(- 


-3) 


1.91(- 


■3) 


(4.19- 


1.52 / tt)(-3) 


6 


0.856(1 


-m)(- 


-2) 


0.43(- 


-3) 


1.07(- 


-3) 


(2.36 - 


0.86^) (-3) 
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TABLE VII: The values of characteristics r^ for £ = and r ni for I > and rt defined by Eqs. (58), (55) and (56) for 
Z = 2 and Z = 10. 





Z = 


2 


Z = 10 


State 




r d 




r d 

' n 


2s 


1.01 




1.00 




2p 


0.68 


1.69 


1.38 


2.38 


3s 


1.28 




1.21 




3p 


0.77 


2.18 


1.68 


2.94 


3d 


0.13 




0.05 




4s 


1.33 




1.31 




4p 


0.80 


2.33 


1.80 


3.18 


4d 


0.17 




0.07 




4/ 


5.0(-3) 




2.0(-3) 




5s 


1.39 




1.33 




5p 


0.82 


2.40 


1.86 


3.27 


5d 


0.19 




0.08 




6s 


1.40 




1.34 




6p 


0.83 


2.43 


1.89 


3.31 


6d 


0.20 




0.08 





